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Abstract 

The applicability of the space-time formulation of the gluonic sector of QCD in 
terms of the Polyakov worldline path integral, via the use of the background field 
gauge fixing method, is extended to multi-gluon loop configurations. Relevant master 
formulas are derived for the computation of effective action terms. 
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Space-time based approaches to QCD, employing first quantization methods, have estab- 
lished a viable alternative for expediting perturbative computations in the theory. The first 
attempts in this direction have used strings as underlying agents [1-3] through which a set of 
prescriptions, known as Bern-Kosower rules, were successfully applied at the one gluon loop 
level. Subsequently, Strassler [4] established the possibility to reproduce the whole approach 
by basing the formahsm on worldline paths instead of strings, see also Refs [5,6]. More 
recently, Sato, Schmidt and collaborators [7,8] have extended the worldline considerations 
to the two loop level. In this letter we apply Polyakov's version [9] of the worldline cast- 
ing of QCD, which we have been consistently pursuing for some time [10-12], to multi-gluon 
loop configurations and formulate specific rules for constructing the relevant expressions that 
represent them. More than just as a strategy pertaining to perturbative calculations this 
approach provides a framework within which a systematic treatment of the propagation of 
particle modes accommodating Wilson lines (loops) can be achieved. 

As in a recent paper [12], in which we presented one gluon loop master formulas for 
effective action terms, we shall formulate the worldline path integral for the gluon by em- 
ploying the background gauge fixing method [13], which is based on a splitting of the gauge 
field potentials by setting = + -B^" and treating the modes as fully fluctuating 
in a background represented by the fields 5^. For the present purposes the latter will be 
considered as classical. The idea, however, that in a quantized context they could be used as 
'carriers' of non-perturbative physics in QCD presents interesting possibilities on which we 
shall devote some comments in the end. In any case, they facilitate the gauge fixing through 
the condition 

(rt^ + gr"'B^^)al ^ Df{B)al = 0. (1) 

Our starting point is the generating functional (partition function in Euclidean space- 
time) given, for the gluonic sector of the theory, by 
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where c, c denote ghost fields, Si is the interaction part of the Yang-Mills action for the a 
fields, Zq is the classical contribution, iG'^J is the gluon propagator and 

R^\a)^{D'r-rgXf'^'< (3) 

Propagators of particle modes (field quanta) according to the Polyakov worldline path 
integral have the following generic form, see Refs [10-12], 

iG^%x,xrL = JdT I Dx{t)exp J dtA (4) 

x{0)=x' ,x(T)=x \ / 

where s = 0, 1/2, 1, ... stands for the spin of the propagating particle entity and m,n are 
'Lorentz generator' indices for each given value of s, e.g., mn ixv for spin 1. Furthermore, 
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$W„[C^^'] is the so-called spin- factor for the particle mode whose generic form is {P stands 
for path ordering) 



^^l[C'^''']^Pexpr- jdtJ-uA ; x{0) ^ X, x{T) ^ x' 



(5) 



with ujpa the tensor ^[xp{t)xcr{t) — Xcr{t)xp{t)] whose indices are contracted with those of the 
generators, namely J^'^ and where [/[C^^ ]"'' is the Wilson line element, specified in terms 
of the background gauge fields, i.e. 



/ T ^ 
C/[C^^']"'' = Pexp ig J dtx{t) ■ B 



ab 



(6) 



In the particular case of the pure gluonic sector, to which we shall exclusively devote our 
attention in this paper, we have {J^vY" (^p/z^o-i/ — ^pi/^a/i)) while for the propagation of 
the ghost modes the spin factor is unity. For completeness, we mention that in the case of 
spin-1/2 modes one has {JmnY'^ l^a^; where the a, P are Dirac spinor indices. 

The above set of formulas suffice to conduct any perturbative, at least, computation 
in the theory. To be concrete, we can produce prescriptions for combining the worldline 
carriers of gluons and ghosts thereby arriving at the worldline version of the Feynman rules 
for the corresponding vertices. Given that the gist of any such computation are the IPI 
configurations, we shall, in the present paper, restrict our considerations to effective action 
functionals, whose one- loop contribution term has been confronted in Ref [12]. For two- 
loops we must separately deal with the configuration having a single 4-gluon vertex and the 
one with two 3-gluon vertices. For the first case we obtain, upon inserting the Polyakov 
worldline path integral expression for the propagators involved and designating by Ci and 
C2 the respective loops that compose the 'figure-8' formation. 
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Xi{0)=Xi{Ti)=x 
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with 



(8) 



the well known 4-gluon vertex. 

In order to deal with the second configuration, which involves a derivative operation at 
each vertex, we first observe that 



D;%zG^^^{x,x')% ^\]dT j Dx{t)x,{T)exp i-\ J dtA 



x{0)=x',x{T)=x 



(9) 
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and similarly, with a relative minus sign, when the covariant derivative acts with respect to 
x' . We now face a situation where the relevant configuration has three, open line, branches. 
The following organization is adopted. Place x andx' on diametrically opposite points on a 
circle. Go, always from x' to a;, in three different ways. Once along one semicrcle, contour 
Ci, once along the diameter, contour C2, and once along the other semicircle, contour C3. 
Following the convention that velocities are positive at the end point we determine, after 
some algebraic manipulations. 



^ 2 3! 
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Dxi{ti)exp 



Xi [Q)=x' ,Xi (Ti )=x 



^^mjcf^^m^^^ (10) 



where the first vertex reads 



+5.4*3(^3) - X2{T2)U, 



(11) 



while the second has a relative minus sign -accounting for opposite velocity orientations- and 
makes the replacement — > 0. 

The last specification needed is for the the two loop configuration which involves vertices 

(2) 

with ghosts. Choosing our paths the same as for , with Ci and C2 traversed by ghosts 
and C3 by gluons and once determining the action of the corresponding covariant derivative 
operators on the ghost field propagator, in analogy to Eq. (9), the following result is obtained 
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Xi{Q)=x' ,Xi{Ti)=x 



x<l>W[Cf">["l[C2^^>[°l[C'r1(V?)?'(K'')5^'^^ (12) 



where 



x\ahc 



(13) 



with 1^2 and Ti ^ for the second vertex. 

One observes that the expressions for the vertices are related to those of the second 
quantized (field theoretical) formulation of QCD via the correspondence p This 
occurence explicitly illustrates the space-time character of our approach to the theory. 

Prom what we have so far presented the following rules follow for the space-time approach 
to perturbation theory: 

•Draw all diagrams (IPl for our present purposes) in terms of worldlines for the particle 
modes which employ the allowed vertices. The analytical expressions for the latter are 
furnished by Eqs. (8), (11) and (13). 

• Por each line segment between adjacent vertices assign the appropriate spin factor and 
Wilson line(loop). 

• Each diagram is assigned a combinatorial factor which coincides with that of the cor- 
responding conventional Peynman diagram. 
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The perturbative computation of the effective action, at two loop level, proceeds once we 
expand the Wilson exponential factors entering the corresponding expressions. Keeping the 
Mth order, in of the expansion means that one determines up to M+2-order contributions 
to the effective action. Before taking this step we find it convenient to recast our two loop 
expressions in a form wherein the various worldlines entering them are described in terms 
of parameters running along the traversed (Wilson) paths as opposed to being characterized 
through the points of reference x, x' . In this way one circumvents zero mode problems due 
to translational invariance while at the same time achieves a form which, in the Feynman di- 
agrammatic language, amounts to directly reaching expressions at the stage where Feynman 
parameters have been introduced. 

The new form of F^^^ becomes, once making the redefinition ^2 — 7\ + ^2 — ^2 and the 
replacements Ti = si, Ti + ^2 = s, 

OO s 

= ds J dsij Dx{t)S[x{s) - x{0)]S[x{s) - x{si)]exp 



where we have used the notation ${i^J,[C**'] = Pexp i ^ J dt J ■ uj\ and similarly for C/[C**']"^. 

Also, we have taken into account the fact that the contour segment C2 of Eq (7) is now being 
reversely traversed since, on account of the uniform parametrization, one branch is given a 
clockwise and the other a counterclockwise orientation. 

For the configurations with the pair of 3-vcrtices the called for changes are Ti+T2+ts 
and t2 ^ Ti + T2 — t2 while the corresponding replacements are Ti = si, Ti + T2 = S2 and 
Ti + T2 + T3 = s. Once taking into account direction reversals since, for the sake of achieving 
a single parametrization, the contour is to be traced in a continuous manner, we obtain 

^ 00 S S2 / 1 * 

Fg^^ = g J ds J ds2 j dsi j Dx{t)5[x{s) — x{si)]5[x{0) — x{s2)]exp J dtx^ 
000 \ 

where 

V,[C,x] ^ /({$W r*-^]Trz.$Wr^°] - $W[C^^^]$W[C^^°]}i:^(.i - 0) 
x[i.(0) + x^{s2 - 0)] + [C-l^WlC-'^] + $W [C^°]K(.i - 0)i:.(0)) (16) 

and similarly for F^^^ with 

V,[C,x] ^ = -/$S[C^^1i^(si - 0)x,(si + 0)). (17) 

We now proceed with the expansion of the Wilson exponential. To Mth order we have 



U[C'''Y^ = 5"^ + ig{t(^T^ j dtixiti) ■ B^^ixih)] + 



+ 



1 



-1 ab 



n ^9{tz. 

.n=M 



1 " 

n /^(in-in-l)i(in)-5^"[2:(i„)]+C»(^^+^). (18) 
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Upon writing 



(2^ 



(19) 



and using the equation of motion D'^F^^^{B) — obeyed by the (classical) background field, 
one determines 



~el{q) = {27rfS{q-pa)el + g{27rr6{q-p,-p,){t1,r2 



b c 
.a \bcc. ^ ' P'^^t. 



(Pb+Pc)^ 



Oig'), (20) 



where the are polarization vectors and the four momenta for the background fields. 
In a generic computation of an n-point Green's function the above expansion serves as the 
means by which a given worldline gluon (muhi)loop configuration is 'punctured' by a network 
of gluonic tree diagrams whose contribution to the overall computation is determined by 
classical field perturbation theory. For the purpose of computing contributions to the effective 
action we are exclusively interested in IPI terms, hence it suffices to set {-B^^[x(t„)]} 
{e^^exp[ipn ■ x{tn)]}. Upon employing the representation [4] 



(21) 



where the CmCn are Grassmann variables, we determine the Mth order term in Eq. (18) to 
take the form 
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+ permutations 


.n=M 




n=M'' i 







In the above expression we have set 



^/i(^n) = Pn,iJ, + ^n^n^^'o^i 



(22) 
(23) 



while the term ^permutations'' refers to all possible rearrangements of the tn- 

As shown in Refs [11,12], via the employment of the Migdal-Polyakov area derivative 
[14,15], for a number of M background gluon field punctures of a given loop configuration 
one obtains the following result for the spin factor 



H'l[C]^H'l[M]=Pexp 



M 



" n=l 

M n-z-l 



M 



n=l 



<7i)ai/ J2 E ^P2a2(^2)0pi<7i(?^l) +■■■ 
n2=l ni=l 



(24) 



with 



.n 



l) = 2Un{s;Pn,. - elPn,,) + Tj^in+lUlUn{el^' - 82+' e;)5{t^+, - t^) (25) 

and where we have set ^m+i = Cm+i = 0. 
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For the generic partitioning, wherein we keep the first N terms (0 < < M) from the 
Wilson loop expansion of the contour C^^^ and the first M — N terms from the contour C^^'^, 
we obtain the following result for the two loop configuration with the four-vertex 



^(2) 

4,(M+2) - g 
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N=0 



-N+l 
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j ds j dsi 


n / dtn 
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n Oitn - tn-l) 



:n=M 



x<l>W [M - iV]$W[A^]Q(^)[{e}, {p}] + permutations, 
where we have denoted (the indication '(1)' signifies the presence of a single vertex) 



(26) 



Q(^) = / Dx{t)d[x{Q) - x{s)]d[x{si) - x{s)]exp J dt^ x^ + ■ x{tn) I • (27) 

V / 

For the three- vertex, two-loop gluon configuration and keeping M — N2 terms from C^^^ , 
N2 — Ni from (7*2*1 g^^^j ^j^g remaining N2 — Ni from (7*^° we get 



^Sgluons 
2,M+2) 
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i?^^-*[{e}, {p}, j]j=o + permutations, (28) 



n ^(^' 

ln=M 



t. 
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N,,N2,-i- 



where the current source j has been introduced for the purpose of facilitating the compu- 
tation of correlators of the type < x^{t)xy{t') > which invariably enter the three- vertex 
computation. Finally,we have set 

r ( 1 } ^ - r 

R^'^\j)^ j Dx{t)5[x{s)-x{Qi)\5[x{si)-x{s2)\exp\--j dt x^ + i^ k{tn) ■ x{tn) + i j dtj-x(tn) 

V "=i 

(29) 

with the '(2)' serving to note that the configuration being considered has two vertex points. 
The ghost-containing configuration results from Eq (28) via the replacement V3 — > V^. 

We are now at a point where we must execute the path integrals. Their gaussian nature 
calls for determining the appropriate Green's functions on each given line contour. It is 
possible to generalize the Green's functions corresponding to propagation of bosonic 'particle' 
modes on unobstructed closed, or open, contours [4-6,12] in (Euchdean) space-time to our 
two- loop situation where one or two intervening vertex points are present. For loops with n 
vertices the Green's function we need can be defined as follows 

(30) 
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with (7(0) (i, t') ^l\t- t'\{s -\t-t'\). 

The corresponding expression for contours with non-coinciding initial and final points 
(open contours), needed for the three- vertex configuration, reads 



with G^°\t,t') = \t-t'\. 

In the above equation Sj^-* is the reduced total time given by 

111 1 

Sl 



(31) 
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and is the reduced partial time 



4") 
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+ ■■■ + 0{tn - S) 



S - Sr, 



(33) 



Using the above expressions one arrives at the following results: 
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1 



(47r)^[s(s - Si)]^/2 



exp 



.m>n 



(34) 



and 



xexp 



(47r)^[s2(s-S2) + si(s2-si)]^/2 



.m>n 



(35) 



For the three vertex computation we need velocity correlators. The result obtains from 
the operation involving the current source j entering Eq (28). One obtains, generically, 

< X,{t)x,{t') > = -5,,dA'&^\t, t') - Y: Htn)Htm)dt[&'\t, 0) - G^'\t, 

n,m 

xdA&'Ht',0)-G^'\t',tJ]. (36) 

Substituting Q^^\ R^^^ and the corresponding correlators in Eqs (26) and (28), respec- 
tively (the latter in the 'ghost' term as well) a set of master expressions is produced which 
call for the execution of integrations over two sets of parameters. We have been able to 
reproduce known results [16] for the M = case and M = 2 cases (forthcoming paper). 

For our closing comments let us begin by saying that the ability to reformulate QCD in 
terms of space-time paths has led to the result of producing master expressions, given in terms 
of two sets of parameters (Grassmann and Feynman) which correspond to multi-gluon loop 
configurations that enter perturbative expansions. Once the overall strategy is established 
the only real determination that needs to be made is the computation of Green's functions 
describing bosonic propagation on paths with a number of 'marked' points satisfying two 
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types of boundary conditions. Even though the examples we have discussed pertain to 
paths that can be continuously traced nothing, in principle, prevents one from treating 
more complex, higher loop, configurations which call for more than a single continuous 
parametrization. The bottom line is that the relevant path integral, whose execution calls 
for determining the Green's functions, is Gaussian, therefore manageable. 

Aside from the execution of multi-gluon loop computations, which call for numerical 
approaches to confront the relevant master expressions as the configurations get increasingly 
complicated, there is another aspect of the Polyakov path integral approach to QCD that we 
find more intriguing. This has to do with the fact that each worldline, in addition to being 
geometrically characterized by the spin of the propagating particle mode through the spin 
factor, is accompanied by a 'cloud' of background gauge potentials on account of the Wilson 
line (loop). Consider a context within which the background fields are quantized in a way 
that they carry non-perturbative physics. This, e.g., can be achieved in a lattice context or, 
alternatively, in lattice inspired continuum schemes, as in Refs [17-19]. One would then be in 
possession of a methodology formulated strictly within the bounds of QCD that will facihtate 
the study of the interplay between the perturbative and the non-perturbative regime of the 
theory. It is to the confrontation of this problem that we intend to direct our future efforts. 
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